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Abstract 

We demonstrate that the modulation instability of the zone boundary mode in a finite (periodic) 
Fcrmi-Pasta-Ulam chain is the necessary but not sufficient condition for the efficient energy transfer 
by localized excitations. This transfer results from the exclusion of complete energy exchange between 
spatially different parts of the chain, and the excitation level corresponding to that turns out to be 
twice more than threshold of zone boundary mode's instability. To obtain this result one needs in far 
going extension of the beating concept to a wide class of finite oscillatory chains. In turn, such an 
extension leads to description of energy exchange and transition to energy localization and transfer 
in terms of 'effective particles' and Limiting Phase Trajectories. The 'effective particles' appear 
naturally when the frequency spectrum crowding ensures the resonance interaction between zone 
boundary and two nearby nonlinear normal modes, but there are no additional resonances. We show 
that the Limiting Phase Trajectories corresponding to the most intensive energy exchange between 
'effective particles' can be considered as an alternative to Nonlinear Normal Modes, which describe 
the stationary process. 

1 Introduction 

It is well known that existence of spatially localized excitations is considerably important in different 
fields of nonlinear physics. Excitations of this type determine elementary mechanisms of many physical 
processes and make a noticeable contribution into thermal capacity [T]. Over the last decades spatially 
localized excitations were studied in the frameworks of two fields of nonlinear dynamics, namely, the 
theory of solitons 0-[5] and the theory of nonlinear normal modes (NNMs) [3 El IB]- The theory of solitons 
was predominantly developed in application to the infinite continuous models described by nonlinear 
partial differential equations. As for NNMs, they can be considered as an extension of linear normal 
modes (LNMs) to discrete nonlinear systems. Development of both theories has revealed the soliton 
excitations in discrete systems (Toda [5] and Ablovitz-Ladik [S] lattices) as well as NNMs in continuum 
models [8 . Existence of the envelope solitons (breathers) allows us to find a common point for both 
theories [5]- 

One of the widely used models of Nonlinear Dynamics is Fermi-Pasta-Ulam (FPU) [10] oscillatory 
chain study of which has become a starting point for discovery of solitons [llj . Long- and short-wavelength 
approximations of the FPU-chain are integrable Korteweg de Vries (KdV) and Nonlinear Schrodinger 
(NSE) equations, respectively [TJ 112) . These equations in the limit of infinite chain have dynamical 
solitons (KdV) and envelope solitons (NSE) as their particular solutions. As for the finite FPU chains, 
their nonlinear dynamics has been intensively studied in terms of NNMs [13] - [23] . 

These studies have shown that for low excitation level a slow energy transfer along such chains is 
realized by spatially periodic waves due to nonlinear resonant interaction between the degenerate NNMs 
corresponding to the same eigenvalue [22 [221 123] ■ The detailed studies of zone-boundary mode (7t-mode) 
in the large enough systems have shown its instability above certain level of excitation [201 [2H [23 [20] • 
Such instability leads to spatial concentration of energy revealed in the form of chaotic breathers, which 
have been studied numerically |20j . However, a simplest model of two weakly coupled nonlinear oscillators 
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(2DoF) allows one to suppose that not one but two dynamical transitions can occur in the systems under 
consideration [57]. In the case of 2DoF system the first of the transitions is caused by instability of 
the out-of phase mode. It is important to note that a complete energy exchange between the oscillators 
can occur both above and below this threshold. But a second threshold exists for a higher excitation 
level and after its crossing the complete energy exchange is forbidden and the energy is confined in the 
initially excited oscillator [57] . The origin of second transition becomes clear if we consider the system's 
phase plane with using appropriate coordinates in which NNMs are depicted as stationary points. The 
loss of stability of out-of-phase NNM leads to creation of two new NNMs corresponding to weak energy 
localization in one of the oscillators. The domains of attraction of these two points are bounded by 
the separatrix passing through the unstable stationary point, and the trajectories providing the energy 
exchange are outside of separatrix. The trajectory corresponding to complete energy exchange is termed 
as Limiting Phase Trajectory (LPT) [27] . It has two branches which limit the domains of attraction 
of unstable (out-of-phase NNM) and stable (in-phase NNM) stationary points, respectively. While the 
excitation energy grows the domain of weak localization expands and the area between separatrix and the 
LPT diminishes up to complete degeneration of latter (when the separatrix coincides with the LPT). At 
this instant the complete energy exchange is forbidden and the energy becomes confined in the initially 
excited oscillator. Let us note that all these results could be naturally formulated in terms of the particles 
(oscillators) rather than of NNMs. 

In the paper we would like to show that similar behavior is common for the nonlinear systems with 
the resonant interaction of NNMs and the analysis based on the LPT concept is a proper tool to clarify 
the fundamental problem of transition from energy exchange to energy localization in a finite (periodic) 
FPU chain. For this, one needs in far going extension of beating concept to the finite oscillatory chains. 
As to best our knowledge such an extension has not been elaborated till now. We have shown that 
introducing the LPT concept leads naturally to the notion of "effective particles" consisting of more than 
one real particles. Then the energy exchange in finite oscillatory chain can be described as beats in the 
system of two weakly coupled "effective particles" similarly to those for two real particles in 2DoF system. 

Our second goal is to fill a gap between the understanding of the origin of the vibration energy 
localization in two-degree-of-freedom and infinite systems. Namely, we derive the conditions of the tran- 
sition to localization in a finite periodic FPU chain with more than two particles. It is worth mentioning 
that the FPU model can be mapped approximately to a DNLS (discrete Ablowitz-Ladik equation) with 
norm conservation ISO] ■ This result is significant for study of strongly discrete high-energy breathers 
because the Ablowitz-Ladik model is integrable one [5]. While dealing with such breathers there is no 
noticeable difference between infinite and finite chain [28] [31] [32] . However, for considered low-energy 
excitations the difference is essential: only for finite chains both energy thresholds studied below exist. It 
is necessary to note also that the mechanism of energy localization considered in the paper differs from 
the scenario 'delocalization-localization' in the harmonic chains with embedded nonlinear oscillator [33j . 



We consider the finite (3-FPU system with periodic boundary conditions defined by the Hamiltonian 



where and p, are the coordinates and the conjugate momenta, respectively, and (qisi + i = c\i, 
Pn+i =Pi)j N is the quantity of the particles. 

The transformation to the normal modes representation (Appendix A) leads to the equations of 
motion in the following form: 



2 The model 




(i) 



l>rrL,n— 1 

Here are the amplitudes of NNMs and the eigenvalues cu^ take the form 




(2) 



UJk = O sin( — ), k 



= 0,...,N-1. 



(3) 
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If considering an even number of particles N, the frequencies tu k are bounded by the uppermost 
value = £1=2. The lowermost eigenvalue cuo = corresponds to motion of the chain as a rigid 

body. All eigenvalues in the interval < cu k < 2 arc doubly degenerate. If the number of the particles 
increases, the frequency gap between the highest-frequency mode and nearby modes quickly decreases 
with increase of N . Taking into account the dependence of the chain properties on the quantity of the 
particles, one can consider the value £ = 1/N as a small parameter of the system. 

To analyze the dynamics of the chain in the context of Eqs ([2]) , we introduce the complex functions 



The complex functions with indices (k, N — k) correspond to the pairs of the degenerate modes 
and as it was mentioned above, their interaction provides the energy transfer by periodic traveling waves 
along the chain. A specific time scale for such energy transfer is proportional e 2 . Both the amount of 
the transferred energy and the rate of the transfer depend on the value of the 'angular momentum' of 
the conjugate normal modes [21] [23] : 



G k - l(M^_ k - V^N-k) - U> k (f, N _ kllk - t k n N _ k ), 

The small value of |G k | corresponds to slow energy transfer, but the large amount of energy 
transferred, and vice versa. The limiting value of G k = corresponds to zero rate of energy exchange 
when the energy transfer is absent. 

The situation changes drastically when the gap between the neighbor eigenvalues diminishes and 
there arise resonance relations between the normal modes corresponding to different eigenvalues. Actually, 
when the number of particles in the chain increases, the frequency gap between the 7t-mode and the nearest 
ones decreases as e 2 . Then a linear combination of the 7t-mode with the nearby conjugate ones may be 
considered as a long- wavelength modulation of 7T-mode with the "wave number" 1 /N . 

In what follows, we consider the coupled dynamics of the 7t-mode (k = N/2) and the nearby 
conjugate modes with indexes (N/2 + 1), (N/2 — 1) and linear frequency cUn/2_i = tuisi/2+i • First of all 
it is necessary to note that the difference of frequencies of the 7t-mode and the nearby ones diminishes 
quickly while the number of particle grows: 

7tf— — 11 7T 2 7T 2 

cun ±1 =Qsin {2 N «Q( 1 -_) = CUN/2 (i-_ e 2) (5) 

Therefore we choose cu N n as a basic frequency and take into account the resonant relationship 
between three highest-frequency modes. The multiple-scale analysis of equations of motion (see Appendix 
B for details) leads to coupled equations (j23|) . ([24]) . 

Keeping in mind the problem of energy exchange and transfer we may restrict ourselves by the 
following equations for three high-frequency modes: 

+ ^[(Ixn/2-i I 2 + IXn/2+1 I 2 )Xn/2 + (x^/2-i +X 2 J/2+1 )Xn /2 ] = 
t d *N/2-. _z^ Xn/2 _ 1 +M[( 2 | Xn/2 |2 + | Xn/2 _ 1 |2 ]xn/2 _ 1 +(2x 2 /2 +x 2 V2+ i)x^ / 2-i] -0 (6) 
■ dxN/2 +1 _^ Xn/2+i +M[(2| Xn/2 | 2 + |x N /2 +1 | 2 )Xn/2 +1 +(2x^/2 +X^ )/ 2- 1 )x^ /2+1 ] =0 

There are three reasons for such restriction: (i) if the low-frequency modes are not exciting from 
very beginning they have no effect on the behavior of high-frequency modes; (ii) the high-frequency modes 
are dominant ones if the initial excitation is spatially localized; (iii) our computer simulation data show 
that the effect of low-frequency constituents is inessential even if their amplitudes are close to those of 
high-frequency modes. 

The eqs © correspond to the Hamiltonian 

H = — 2 2~(IXN/2+l I 2 + IXN/2-1 I 2 ) + T^L2;(IXN/2+1 I 2 + IXN/2-1 I 2 ) 2 + 2"(Xn/2+1Xn/2-1 _ Xn /2+1 Xn /2-1 ) 2 
(XN/2+1Xn/2 +Xn/2+iXN/2) 2 + (XN/2-1Xn/2 + Xn/2-1 XN/2 ) 2 ] 
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(7) 

Eqs ^ possess an additional integral of motion, which is said to be the total 'occupation number': 



X = Ixn/21 2 + IXn/2-i I 2 + IXn/2+i I 2 = const (8) 
(we omit the index of X because in following we will consider the highest-frequency modes only) . 

3 The 'effective oscillators' description 

To analyze the system, we turn from the 'wave representation' to the 'particle representation'. Really, 
each of the modes considered represents a delocalized solution with the spatially uniform (7t-mode) or 
periodic (conjugate modes) energy distribution. However, as it was mentioned above, a combination of 
the modes gives raise to a spatially non-uniform solution, corresponding to the energy localization in the 
one half of the chain. So, choosing a linear combination of the considered modes we make obvious an 
analogy with the system of two nonlinear oscillators, where each half of the chain is considered as an 
'effective oscillator'. This implies the transition from the 'wave' to 'particle' representation. In such a 
case the beating in the system of the 'effective oscillators' manifests as the process of periodic energy 
transfer from one half of the chain to another one and inversely with time scale z 2 . To this end, we 
introduce the 'particle' representation through the variables ibi , i|>2, <P, defined by the following formulas 



^1 = 7f (XN/2 - (Vl -C 2 XN/2-l + CXN/2+1 )) 

^2 = 75 (XN/2 + (Vl -C 2 XN/2-l + CXN/2+1 )) , (9) 
tP = (cXN/2-1 - Vl -C 2 XN/2+1 ) 

where c is a constant defined by initial conditions (0 < c < 1 ) . Such transformation preserves the 
occupation number in the form X = |ibi | 2 + ji|>2l 2 + M 2 - 

We consider first a particular solution corresponding to <p = 0. Since the occupation number X is 
the integral of motion, the variables ipi and t\>2 can be expressed through the 'angular variables' [M] : 

xbi = VXcos9e i6l ; ib 2 = \/Xsin9e iS2 . (10) 
Then the Hamilton function Q takes the following form in the terms of 'angular variables': 



H(9,A) = ^[27(3X-16tt 2 +2(87T 2 -3(3X) 
x cos A sin 20 — 3|3X(8 — cos 2 A) sin 2 29], 



(11) 



where A = 62 — 61 . 

The analysis of phase-plane portrait of the Hamiltonian (II 1 D shows that there exist two stationary 
points corresponding to the pure 7t-mode (9 = 7t/4, A = 0) and the "pure" mixture of conjugate modes 
(9 = 7t/4, A = 7t), respectively, at a small value of the occupation number X (Fig. la). If we consider the 
phase trajectory associated with the transition from the state i|>i (9 = 0) to the state i|>2 (9 = 7t/2), then 
we find that it corresponds to beating in the system of the "effective oscillators" . Since this trajectory 
represents an outer boundary for a set of trajectories encircling the basic stationary points (Fig. la), we 
refer to it as the Limiting Phase Trajectory (LPT) [27] . The motion along this trajectory leads to the 
complete reversible energy transfer between two "effective oscillators" . 

If the excitation level grows and attains a critical level the 7t-mode becomes unstable that leads to 
the appearance of two new stationary points (Fig. l.b) with the folliwing coordinates : 

(a) A = 0, 9 = larcsin(^M) 

(b) A = 0, 9 = f -^ajcsin( 8 * 2 - 3 x px ) ' 

The threshold of instability of the 7t-mode is equal to X c = 7t 2 /3(3, which is in a good accordance 
with the estimation obtained in the framework of the "narrow packet" approximation |251 120j . 
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Figure 1: (color online) The truncated phase-plane portrait before (a) and after (b) threshold of insta- 
bility of the p-mode; the (red) dashed and (blue) long-dashed curves correspond to the LPT encircling 
7t-mode and the separatrix, respectively. 




=JT — 71/2 ^11 A n 



Figure 2: (color online) Phase-plane portrait after coincidence of the LPT and the separatrix; (black) 
long-dashed lines, and (blue) short-dashed lines correspond to the separatrix, passing through the unstable 
stationary point, and the transient trajectories, respectively; the (red) dotted line depicts the boundaries 
between the families of transient trajectories and closed ones encircling the stationary points. 

Any trajectories lying in the neighborhood of stationary points (11 2[) correspond to weakly localized 
solutions, for which the energy of one half of the chain ("effective oscillator") only slightly exceeds the 
energy of the second one. At the same time, the trajectories starting beyond this vicinity and passing 
through the points ibi (9 = 0) or ibi (9 = 7t/2) correspond to the combinations of modes with ap- 
proximately equal energies that means that an initial excitation of one "effective oscillator" entails the 
complete energy exchange between both ones, i.e. the transition from the state ibi to the state ib2 and 
inversely. This implies that a possibility of complete energy exchange between different parts of the chain 
exists for the excitation level exceeding the instability threshold X c . 

A growth of the amplitude X entails an enlargement of the domain encircled by the separatrix 
passing through an unstable stationary point (out-of-phase NNM); at last the separatrix coincides with 
the LPT. At this instant topology of the phase plane changes drastically (Fig. 2). Namely, any energy 
exchange between the "mixed" states ibi and il>2 disappears; this implies that a trajectory starting at 
a point corresponding to 9 < 7t/4 (or 9 > 7t/4) and for any A cannot reach a point corresponding to 
9 > 7t/4 (9 < 7t/4) (excepting all the trajectories surrounding the in-phase stationary points A = ±7t and 
bounded by the separatrix crossing the unstable point (9 < 7t/4, A = 0)). Therefore, the energy initially 
concentrated near the states ibi or i|>2 remains confined in the excited "effective oscillator" (Fig. 3). 

An energy threshold associated with the above mentioned localization can be found from the con- 
dition of equality of the energy corresponding to the LPT and the energy at the unstable 7t-mode. It is 
seen in Fig. 1 that the LPT goes through the points 9 = n/2 and 9=0. This means that 

H(9,A)| LPT = ^(27PX-I67t 2 ) 

Since the energy of the 7t-mode is equal to zero, it is easy to calculate the respective occupation 
number X\. oc = 167t2/27(3 and the energy of the chain Xi oc — 167t2/27|31M. It now follows that above the 
excitation level Ei oc we can observe the localized vibration excitation (a breather). 
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Figure 3: (color online) Transition "beating - localization" in the FPU-chain with 10 particles. The 
initial conditions correspond to the state iJji , and the levels of excitation are X < Xi oc (a) and X > Xi D c 
(b). E and n - the energy (in the relative units) and the number of particles, respectively, t - the time 
(measured in the periods of 7t-mode - To = 27t/cu N /2)- 



G 



If the initial energy is concentrated at the state ibi or il>2, the representing point in the phase 
plane moves along a trajectory encircling the respective stationary point. Then the temporal evolution 
of the breather corresponds to the regular variation of its profile (the "breathing" mode of the localized 
excitation). 

The period of breathing can be calculated as 



f 1 dA 

T = odt= -5-0—7- — , 
J e 2 J dA/dx 2 

where the last integral is taken along the LPT. 

Contrary to the "breathing" breathers corresponding to the motion along the LPT, the stationary 
points (|12l) determine new normal modes with an invariable non-homogeneous energy distribution, or the 
breather-like excitations. Although these modes exist at any X > X c the true threshold of localization is 
equal to X\ oc because of a possibility of the complete energy exchange is preserved for motion along the 
LPT till X = X loc . 

Thus we get the solution (|12p corresponding to immobile breather. Finally we want to clarify the 
nature of traveling breathers. We recall that the above results have been obtained under assumption 
(p = that corresponds to equal amplitudes of the conjugate modes. Now we assume that the amplitudes 
of the conjugate modes are slightly different, that is we consider small enough cp ^ 0. One can show 
that the equations for functions ipi and i\>2 include only quadratic terms depending on <p. Therefore, in 
the framework of the linear approximation, if the value cp 7^ is small enough, there is no qualitatively 
change in the phase-plane portrait in figures (1, 2). The behavior of cp in the vicinity of any of stationary 
points (fT2"| is described by following equation: 

. dep _ tt 2 + 3 |3 X[2sm20(p + (3 + sm2 ee 2iS| )<p*] = 0. 
dT2 2 1 6 

The respective eigenvalue 



v/891 |3 2 X 2 + 120tt 2 |3X - 384tt 4 
A = ± 56 

is imaginary if X < Xi oc , i.e. if the amplitude of excitation is less than the threshold corresponding 
to the coincidence of LPT and the separatrix. 

The distribution of the vibration amplitude along the chain can be characterized by the mean 
squared particle displacements < q 2 >: 



2 <^N/2 



■2n/w N/2 

q 2 (t,T 2 )dt. 





For cp = 0, then there exists a single maximum of the vibration amplitude at j — N/2. It is easy 
to demonstrate that, if the eigenvalue A is imaginary, then a small nonzero cp in the equations of motion 
leads to small oscillations of the breather center. At the same time, appearance of a real part in the 
eigenvalue A for X > Xi oc leads to the directional motion of the breather (Fig. 4). 



4 Discussion 

Concluding, one can state that the energy localization in the finite nonlinear system like the FPU lattice is 
caused by the resonance interaction of nonlinear normal modes which are close to the spectrum boundary. 
At that, the adequate description of the process has achieved in the terms of "effective oscillators" rather 
than in terms of NNMs. There are three levels of the system excitation, at which the behavior of the 
system is obviously different. At a low excitation energy the resonance interaction of modes leads to a 
phenomenon that is similar to beat oscillations in a system of two weakly coupled nonlinear oscillators, 
i.e. to periodic energy exchange between two "effective oscillators". Both the complete and partial 
energy exchange can be realized. The loss of stability of the zone-boundary mode leads to weak energy 
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Figure 4: (color online) Energy map for static (X < Xi oc ) (a) and traveling (X > Xi oc ) (b) breathers. 

localization, when only a small excess of the energy is confined in the one "effective oscillator" . Such 
energy excess increases with the growth of the excitation level. As this take place, the possibility of 
complete energy exchange is preserved in the attraction domain of "out-of-phase" oscillations, while 
both the complete and partial energy exchange can occur in the attraction domain of the "in-phase" 
oscillations. The center of the energy localization can oscillate under the influence of small perturbations. 
When the excitation level reaches the second threshold correlated with the degeneration of the area 
bounded by the LPT the beating converts into strong localization. This leads to the energy confinement 
in only one "effective oscillators" (i.e. in one part of the chain). It is important that this transition 
involves both "out-of-phase" and "in-phase" trajectories, i.e. the localization does not depend on the 
initial "phase shift" A. A weak localization and non-complete energy exchange are still possible but their 
realization requires some special initial conditions. One of the key changes after the second transition is 
the transformation of small oscillations of the center of the localization area into a directed motion. Also, 
since the domain of strong localization does not correspond to any stationary point in the phase space of 
the system, the profile of the energy distribution along the chain demonstrates the periodic variations. 
A further growth of the excitation does not result in any qualitatively changes in the systems behavior, 
which might increase or decrease the areas corresponding to both partial energy exchange and weak 
localization. It should be note that increase of the particles number in the chain involves the new modes 
into resonance interaction with zone boundary mode because of spectrum crowding in the high-frequency 
range. It does not lead to any noticeable changes except some "narrowing" of "effective oscillator" . In the 
limiting case of infinite chain the "effective oscillator" transforms into well-known breather. At that, both 
energy thresholds discussed above tend to zero. It means that the resonant nature of energy localization 
can be better clarified while dealing with the finite chain. 



5 Conclusion 

Finally, we have revealed the origin of localization occurring due to resonance interaction in the range 
of the spectrum crowding in relatively small FPU chains. Using the LPT concept and the notion of 
"effective particles" has allowed us to clarify the origin of transition from energy exchange to energy 
localization and transfer in finite periodic oscillatory chains. We highlight the role of the size of systems 
under considerations. Since the gap between high-frequency modes decreases as 1 /N 2 the resonant re- 
lations are already realized even if the chain contains a small enough number of particles. Actually, the 
phenomenon of strong energy localization is seen clearly in the system with 8-10 particles, and the com- 
puter simulation data are in the excellent agreement with the analytic results. We would like to underline 
that the number of resonance relations increases and both discussed thresholds tend to zero when size 
of the chain grows. Therefore the origin of energy localization can be effectively treated if considering a 
finite chain with small enough number of resonances. Two non-zero thresholds exist for any finite chain 
and the only quantitative change in the case of large number of particles is decrease of the "size" of 
the "effective particle" because of the influence of new resonances. Small perturbations of the "effective 
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particles" behavior owing to the modes, which can be ignored in the resonance approximation, become 
visible on very large time intervals (several thousands of the oscillation periods). We have observed that 
a long-run temporal evolution leads to the formation of localized excitations with the behavior similar to 
that of chaotic breathers. 
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Appendix 
A The normal modes 

The finite |3-FPU system with periodic boundary conditions described by the Hamiltonian ([1]) can be 
written in the terms of normal modes which are defined by the following linear canonical transformation 

N-1 

qj = Y_ (13) 

k=0 

where 



o- j)k = J=[ s hi(^) + cos(^)], k = 0, 1, N - 1 . 



Transformation (TT5)) allows us to present the quadratic part of the Hamiltonian as the total energy 
of the independent oscillators: 



N-1 



H 2 = X 



k=1 



(14) 



Here £,k and r|k are the amplitudes and the momenta of NNMs; the coordinate f,o, associated with 
motion of the center of mass, is removed from (TT4"]) . The eigenvalues u>k are defined by eq. 

If the number of particles N is even one, the frequencies a>k are bounded by the highest-frequency 
cujsi/2 = O = 2. The zero eigenvalue a>o — corresponds to motion of the chain as a rigid body. All 
eigenvalues in the interval < o>k < 2 are doubly degenerate. The density of eigenvalues a)k is the 
essentially nonuniform function of k: the lowest values (small k) are separated by large gaps while the 
distances between highest frequencies (k near the N/2) are small enough and they diminish quickly with 
the growth of number of particles. 

The 4th-order part of the Hamiltonian (TTJ) for the periodic (3-FPU chain has the following form [T3] : 



H 4 = Y_ W k <^l<^mUJ n Ck ) l ) m ) n£,k£,l£,m£,n, ( 15 ) 

k,l,m,n— 1 

where 



Ck,l,m,n — "Ak+l+m+n + Ak+l-m-n 

(16) 



"Ak-l+m-n + Ak-l-m+n 



A r = 



(-1) T , if r = mN, m e Z 
otherwise. 



It now follows from fH|) . (|15J) that the equations of motion can be written in form $Z§. 
To analyze the dynamics of the chain we introduce the complex functions 
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= — (uk + io) k £,k), = — (rik -icu k £, k ) 



1 

In terms of these variables Eqs ([5]) can be rewritten as 



B The multiple scale expansion 

To consider the processes on the time intervals greatly exceeding ^/wy_ we employ the multiple scale 
procedure. We introduce the time scales: 

T =t, Tt^ETo, t 2 = £ 2 t ,..., (18) 

where the "fast" time To corresponds to the original time scale of the system, while the slow times 
Ti , T2, etc. correspond to the slowly varying envelopes. 

We construct an asymptotic representation of function Wy_ in the form 

^k = Ve(<Pk,i + £(Pk,2 + £ 2 <Pk,3 + •••)• (19) 

Taking into account the relations (fTS)) and (|19D we get from Eqs (|17[) in the lowest-order by small 
parameter e: 



+ w n cp n 1 = 



^ .IT?.,: ' ■ (20) 



i ' g'' 1 +o)n_ 1 ( P n ±11 =0 
Taking the solutions of Eqs ([20)1 in the form 



CPN/2,1 =Xf,ie iaTo 
cpN ±1)1 =XN ±1)1 e iQT ° 



we get the discrepancy in the second of eq. (f20|) that is the order of t 2 and it has to be taken into 
account in the approximation of next orders. The amplitudes Xnj i X-M-±i i do not depend on the 'fast' 
time To. 



3<P N 2 9<p N , 

' i I 1 h CU n (D n 7 + i 5 J' = 

.3/2 . J 3t ~^~' z 9t i 

— ^— + (0 N _! (p n ±1 >2 + 1— g^— = 



Taking 



cpN ±1)2 =XN ±1)2 e^° ' 



we get: 



3xn , 
i — ^7—^ = 

9x n . , , 

i J* ■ = 

OTi 
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The difference between ix>n_ and oin_ ] is negligibly small in this order by the small parameter e. 



( 3( p j 



,V2 



3cp j 



,2 



3(p j 



Tll!m!n=l C ^ ) l,m ) n[ ( Pl,1 < Pm ) lCp 1 x,1 ~ 3<Pi,l <p m ,l <Pn,1 + 3t Pl,1 <Pm,1 <Pn,1 ~ <Pl,1 <Pm,1 <Pn,l] 

3 f N ±1 2 3<P N ±1 3 _2 

+ t J hCUN._ 1 (pN. ±1 3 ^Hcpj 



■ 9<Pn/2±i , l , ■ 

1- "T X 



k ^rLlIm!n=l C ^-±l,l,m,n[«Pl,1 < Pm l 1<Pn,1 ~ 3(pl,1 <Pm,l <Pn,1 + 3 <Pl,1 <Pm,1 <Pn,1 ~ «P*,1 <Pm,1 <Pn.ll 



Taking 



cpN.,3 =XN 3 e 



[ qjN ±13 =XN ±1)3 e 
we get: 



J30 r-N-1 



Vll,m,n=l C f ,1 ,m,n[<PM <Pm,1 <Pn,1 - 3<pl,1 <p m ,1 <p*j + 3cp t)1 (p^, <p* (1 - (p^ Cp^-, (p'^Jc" 
9 X. 



L ±1<1 . 3 Xn ±1i2 2 
1 »~ + X ^xl T%±1,1 



3T2 



f ljm!n=1 C^-il.l.m.nlVl.KPm.KPn.l ~ 3(pl,1 <pm,1 <Pn,1 + 3( Pl.l <Pm,1 <Pn,1 ~ <Pl,1 <Pm,l <Pn,lle 



(21) 



To remove the second-order constituents from the last equations one should integrate them with 
respect to the times To and Ti . To do it we have to know the manner in which the variables <pk depend 
on the 'fast' time. We can perform the similar procedure for all variables cpk- 



J/2 . 



9t 



•U>lc<Pk,l =0 



I 9t + t^icCpN-lc,! — 



The obvious solutions of the last equations are: 



<Pk,1 = Xkje LLi 
<PN-kJ = XN-k,1 e 



lO) k To > 



where the amplitudes Xkj and Xn— k,l do not depend on the time To. In the next order of £ we get 
the following relationships: 



,3/2 . 



i a i£ £_^ + a , k(pN _ k)2 + i^-^ = 



<Pk,2 = Xk,2 e 
(PN-k,2 = XN-k,2S 



lO) k T > 



and 



3Xk 



9ti 

9XN-k,i = 
3t, 



The amplitudes Xk,i , XN-kj , and Xk,2, and XN-k,2 do not depend on the times To and Ti . Now we can 
integrate the equations (j!?Tj) with respect to the time Tq. The requirement of absence of secular terms 
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means that only the terms with zero- value argument of exponent make a contribution in the sum in Eqs 
(f2"Tj) . Finally, after the integration of obtained equations with respect to the time Ti we get: 



9 X N 1 TO XI 1 

i ^T + T Q Il,m,n=l CN ilm)n Xl,lXm,lX;, 1 6(cWl + CU m -U; n -n)=0 

9 X N ± , , 2 30 N— 1 

1 ?t 2 ' ~ T- Q XN/2±1,1 + T 2-l,m,TT=1 C ^±1 ,l,m,nXl,lXm,lX^l + ~ ^ ~ °J = 

It is necessary to note that the frequency spectrum of a finite FPU-system is essentially discrete one 
and no exact resonances exist except those for the degenerate modes with indices k, N — k. Therefore 
only approximate resonant conditions can be satisfied for the highest-frequency modes, as it has been 
pointed in the Introduction. Taking into account these rules we find the final form of Eqs (|22[) (we omit 
the second index in the subscription of main-order amplitudes Xk): 

i w + ^KLlLl 1 bal 2 + Ixn +1 1 2 + Ixn | 2 ) X n + (xi +1 +x|_! )X*n] = o 
^r- - z 4 q xn +1 + ^{2 Lit. 1 IXil 2 X f + i + K2|X^I 2 + Ixn +1 i 2 + 

Ixn-iI 2 -|xiI 2 )Xn+i +(2 X 2 n +X 2 n_Jx*n +1 -(xiXn_i + Xn-iXi )x« -,]} = (23) 



i^T 1 - t 1 ^-, + ^{2 Lit. 1 IXil 2 X f -i + [(2| X n I 2 + | X n_, I 2 - 
Ixn-iI 2 + IxiI 2 )Xn_ 1 +(2x 2 n +X 2 n +1 )x*n_! - (XiXn-i + Xn-iXi )Xn +1 ]} = 0. 



It is easy to show that the dynamics of the low-frequency modes xi and Xn— i is defined by the 
next equations: 



+ ^[(Llti Ixil 2 + Ixi l 2 )xi +xi^x] + (Ixn/2+i I 2 - Ixn/2-1 l 2 )xi- 

(Xn/2-i Xn/2+i +Xn/2+iXn/2-i )Xn-i1 =0 
(Xn/2-i Xn/2+i + Xn/2+i Xn/2-i )Xil =0 



(24) 



First of all we can see that Eqs ([231) , (|2"4")l couple the high- and low- frequency modes with the pairs 
of indices (k, N — k) and (N/2 — k, N/2 + k). This result is in accordance with the resonant normal form 
of FPU chain with even number of the particles [T7]. Moreover one can find that Eqs (|23p . (|2"4"]l posses 
the integrals of motion which are an additional ones to the integrals of energy end momentum: 

X 4 = Ix#| 2 + Ixn +1 | 2 + Ixn_ 1 | 2 



Xi 2 =IxiI* + Ixn-iI 2 



The existence of these integrals allows us to conclude that the sum over all the modes is a constant: 
Li^i 1 Ixil 2 — LwfXi = const. The presence of such terms in Eqs (|2"4"]l leads to the frequency shift 
5cuk = (3|3u>k/4) Lil^T 2 X\- In such a case we can redefine the amplitudes: Xk — > Xk6 l6tUkT2 , and extract 
the respective terms from the equations. 
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